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Abstract. We give explicit numerical values with 100 decimal digits for the Mertens con- 
stant involved in the asymptotic formula for Y, ^/p ^nd, as a by-product, for the Meissel- 

p=a mod g 

Mertens constant defined as Y^p=a mod q 

(log(l -l/p) + l /p), for ^ e {3, . . . , 100} and {q,a) = 1. 
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1. Introduction 

In this paper we use the technique developed in [|71 to compute the constants M{q, a) involved 
in the following asymptotic formula 

i^loglogx / 1 X 

p=a mod q 

where x — > +o°, and the so-called Meissel-Mertens constant 

B{q,a):= £ (log(l-i) + i 

p=a mod q P P 

where, here and throughout the present paper, q >3 and a are fixed integers with (q,a) = I, p 
denotes a prime number, and (p{q) is the usual Euler totient function. In fact we will see how 
to compute M{q, a) with a precision of 100 decimal digits and we will use the results in [7J to 
obtain the values for B{q, a) . 

To do so we recall that the constant C{q, a) studied in |[6l|7l is defined implicitly by 

:^ n ('-^)-(i|fi|iT('+''(')) 

p=a mod q 

asx—^ In O we proved that 

1 ^ a{p;q,a) 

P' 



c(^,a)*p(^)=^-^n(i-- 



where a{p;q,a) = (p{q) — I if p = a mod q and a{p;q,a) = — I otherwise, and yis the Euler 
constant. This enabled us to compute their values with 100 decimal digits in fTl . 
Taking the logarithm of both sides in Q we get that 



I .os(l-i)^.ogC(,.„)-».<„.) 



p<x 
p=a mod q 



as X ^ and hence, adding ([T]), we obtain 



M{q,a) = B{q,a) —logC{q,a). (3) 
1 



By ^ and using the results in [|7]| together with the computation on M{q, a) we will explain, 
we can compute the corresponding values for B{q,a) in the same range (and with the same 
precision) for any ^ G {3, . . . , 100} and {q, a) = I. 

We recall that Finch 0\ has computed M{q, a) and C{q, a) in the case q E {3,4} and {q, a) = 

1. 

Acknowledgments. We would like to thank Robert Baillie [[Q who has driven our attention 
to the problem of computing M{q, a) . 



2. Theoretical framework 

From now on we will let % be a Dirichlet character mod^. By the orthogonality of Dirichlet 
characters, a direct computation and Theorem 428 of Hardy-Wright L5J show that 

cp(^)M(^,a)=Y+5-£i+ £ (4) 

p\q ^ X mod q P P 

Xt^Xo 

where 

_ / 1 1 \ 

(5) 



5:=i:(log(l--) + - 



is the Meissel-Mertens constant. Moreover, using the Taylor expansion of log(l —x) and again 
by orthogonality, it is clear that 

^MBM = - E xwE^E^ = - E ME^E^+«(«). (6) 

Xmodq m>2'" p X mod ^ m>2 p t' 

Xt^Xo 

where B{cj), defined as 

B{q):=-T- E — : 

m>2^'^^ {p,q) = \P 

represents the contribution of the principal character %o mod q and it is equal to 

(i„g(,_i,+i)=s-E(iog(i-i)+i; 

(p,g) = l r f r r 

where B is defined in ([5]). Recalling from section 2 of [|7l| that 

cp(,)logC(,,a) = -Y+log-^- £ ml^lx'^ il) 

Xt^Xo 

and comparing the right hand sides of ([4]), ^ and ([7]), it is clear that it is much easier to compute 
M{q,a) than both C{q,a) and B{q,a) since in (HI) no prime powers are involved. Moreover, by 
([3]), we can obtain B{q, a) using M{q, a) and C{q, a) . 

Since in fJl we already computed several values of C(q,a), it is now sufficient to compute 
M{q,a) for the corresponding pairs q,a. 

To accelerate the convergence of the inner sums in (H)), ^ and (|7]), we will consider, as we 
did in [7], the "tail" of a suitable Euler product. Letting A be a fixed positive integer, we denote 
the tail of the Euler product of a Dirichlet L-function as 



p>Aq ^ P ^ 



1 



2 



where x¥'Xo mod q and > 1. Now we prove that 

1 2cM = £^log(L^,(x^M), (8) 

p>Aq P k>\ 

for every integer m > 1. We use the Mobius inversion with a little care, since the series for 
LAqix, 1) is not absolutely convergent. The Taylor expansion for log(l —x) implies that 

k>2 p>Aqk>2n>l "'^/^ p>Aqe>2 k>2 

k\i 

p>Aql>2 ^P p>Aq ^ 

since A'(^) = for £ > 2, and this proves ([8]) for every m > 1. 
Inserting now ([8]), with m = 1, in (01), we have 

(p(^)M(^,a)=Y+5-£i+ £ %W L ^+ I XHi:^log(W%''^)) 

^ X mod g p<Aq P xmoAq k>\ 
Xt^Xo Xt^Xo 

= (p(^) £ -+M(^)+ £ x(a)£^log(L^,(%^fc)), (9) 

p<A^ ^ X mod ^ /:> 1 

p=a mod q x/Xo 

where 

M(^):=Y+5-£-- £ -. 

(p,9)=l 

For A > 1, it is clear that the two sums at the right hand side of the previous equation collapse to 
Y.p<Aq 1 but in (flOl) we will explicitly need the value of the summation over p \ q and hence, 
to avoid double computations, we will use the definition of M{q) as previously stated. 

For C{q,a) the analogue of ^ is eq. (5) of OJ while for B{q,a) it can be obtained arguing 
in a similar way. 

Notice that the Riemann zeta function is never computed at 5 = 1 in ([91), since for ^ = 1 we 
have %^ = X = Xo- To compute the summation over % in dH) we follow the line of section 2 of 
0. 

This means that to evaluate ^ using a computer program we have to truncate the sum over 
k and to estimate the error we are introducing. Let ^ > 1 be an integer. We get 

<?iq)Miq,a) = <?{q) £ -+M{q)+ £ £ ^log(LA,(x',^)) 

p<Aq P Xmodg 1<A:<A: 



p=a mod q Xv^Xo 

k 



+ I mi,^iog{LAq{x\k)) 

X mod q k>K 
X/Xo 

= M{q,a,A,K) +Ei{q,a,A,K), 

say. We remark that B, defined as in ([5]), can be easily computed up to 1000 correct digits in 
few seconds by adapting ^ to the case in which the sum in the left hand side runs over the 
complete set of primes. We recall that Moree [ISJ], see also the appendix by Niklasch, computed 

3 



B and many other number theoretic constants with a nice precision, see also Gourdon-Sebah's 
[HI website. Using the Lemma in [7J and the trivial bound for %, it is easy to see that 

2(A^)i-X^)-l) 

We take this occasion to correct a typo in |[71 in which, in the inequality for E\{q,a,A,K) at 
page 319 there, the factor 2K at the denominator should be read as K^. 

In order to ensure that M{q,a,A,K) is a good approximation of M{q,a) it is sufficient that 
Aq and K are sufficiently large. Setting Aq = 9600 and K = 26 yields the desired 100 correct 
decimal digits. 

Now we have to consider the error we are introducing during the evaluation of the Dirichlet 
L-functions that appear in M(^, a,A,K) . This can be done exactly as in section 3 of [VJ replacing 
km there by k. Let T be an even integer and be a multiple of q. For x¥'Xo mod q and k> I, 
we use the Euler-MacLaurin formula in the following form 

T (.k ^^ V 1 V i-'^y-'Bjix')k{k+\)---ik + j-2) 

where denotes the x-BernouUi number which is defined by means of the n-th Bernoulli 

polynomial 5„(.)c) (see Cohen [2], Definition 9.1.1), as follows 

in which / is the conductor of %. 

Hence the error term in evaluating the tail of the Dirichlet L-functions L^^(x^,^) is 

|.,„,.,^,^,r)|<(|f)-l^ £ >^i^±Iz^^.-- 

^{<^{q)-l)q^Br y ^k+,y..^k + T-2)N'-'-^ 
U{q,K,NJ)T\,^^^^ ^ ^ > 

^ {^{q)-l){K + T-2Y-^q'BT y , 
U{q,K,N,T)NT-iT\ ^J^^ 

^ 2{^>iq)-l)iK + T-2Y-^q^BT 
- {N-\)U{q,K,N,T)N^-^Tl ' 
where Bx is the T-th Bernoulli number and 



U{q,K,NJ):= min min |Lr,^(%^fc)|. 

X mod ^ 1 <A:<A: 
Xt^Xo 



Collecting the previous estimates, we have that 



M{q,a) 



M{q,a,A,K) \E{q,a,A,K,N,T) 



< 



(?{q) (p{q) 

where E{q,a,A,K,N,T) denotes Ei{q,a,A,K) + E2{q,a,K,N,T). 

Practical experimentations for q E {3, . . . , 100} suggested us to use different ranges for A'^ 
and T to reach a precision of at least 100 decimal digits in a reasonable amount of time. Using 
Aq = 9600, K = 26 and recalling that q \ N and T is even, our choice is = ( [8400 /q\+l)q and 
r = 58 if ^ G {3, . . . , 10}, while for q G {90, . . . , 100} we have to use = ( [27720/^J + 1)^ 
and T = 88. Intermediate ranges are used for the remaining integers q. 
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The programs we used to compute the Dirichlet characters mod^ and the values of M{q,a) 
for ^ e {3, . . . , 100}, 1 < a < ^, a) = 1, were written using the GP scripting language of 
PARI/GP f9l; the C program was obtained from the GP one using the gp2c tool. The actual 
computations were performed using a double quad-core LinuX pc for a total amount of com- 
puting time of about 4 hours and 4 minutes. 

A tiny part of the final results is collected in the tables [Ill6] listed at the bottom of this paper. 
The complete set of results can be downloaded from http : / / www . math . unipd . it/~languasc/ Mertens-compu 
together with the source program in GP and the results of the verifications of the identities (flOl) 
and (fTTl) which are described in the section below. 

Moreover, at the same web address, you will also find the values of B{q,a) computed via 
([3]) using the previous results on M{q, a) and the ones for C{q, a) in Q. The use of ^ implies 
some sort of "error propagation". To avoid this phenomenon we recomputed some values 
of C{q,a). A complete report of this recomputation step can be found at the web address 
previously mentioned. 

Moreover, to be safer, we also directly computed B{q,a) using ([61) for ^ G {3,..., 100}, 
\ < a < q and {q,a) = 1. The needed computation time was about 3 days, 6 hours and a 
quarter. By comparing the values of B{q,a) obtained using these two different methods, we 
can say that the values of B{q, a) we computed are correct up to 100 decimal digits. 

Finally, we also wrote a program to compute B{q,a), C{q,a) and M{q,a) with at least 20 
correct decimal digits. Comparing with Q, the main parameters can be chosen now in a 
much smaller way and so we were able to compute all these constants for every 3 < q < 300, 
I < a < q, {q,a) = I. In particular, the needed time on a double quad-core LinuX pc for the 
range ^ e {3, . . . , 200} was about 5 hours and 5 minutes while, for the range q E {201 , . . . , 300}, 
it was about 18 hours. In this case we directly computed B{q,a), C{q,a) and M(q,a) and we 
used ([3]) as a consistency check. 

The whole set of these results can be downloaded at the web address previously mentioned. 



3. Verification of consistency 

The set of constants M{q,a) satisfies many identities, and we checked our results verifying 
that these identities hold within a very small error. The basic identities that we exploited are 
two: the first one is 

£ M(^,a)=y+5-£-. (10) 

amodq p\<l ^ 

This can be verified by a direct computation, taking into account the fact that primes dividing 
^ do not occur in any sum of the type 1^ p<x ^. 

p=a mod q 

The other identity is valid whenever we take two moduli qi and qj with qi\q2 and (a, <7i ) = 1 . 
In this case we have 

n—l 

M{qi,a)= £ M{q2,a + iqi)+ £ - (11) 

J=0 p\qi P 

{a+jqi ,q2) = l p=a mod qi 

where n = q2/q\- 

Equation (fTTl) holds also for B{q,a) with the only remark that in the final summation the 
summand l/p should be replaced by log(l — l/p) + l/p). Concerning (fTOl) . this holds for 
B{q,a) too if we replace l-"Lp\q 1/P with -Lp|,/log(l -\/p) + \/p)). 

The proof of ([TT]) depends on the fact that the residue class a mod q\ is the union of the 
classes a + jq\ mod qj, for j E {0, . . . , n — 1} . li q\ and q2 have the same set of prime factors 
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the condition (a + jqi^qi) = 1 is automatically satisfied, since (a,^i) = 1 by our hypothesis. 
On the other hand, if q2 has a prime factor p that qi lacks, then there are values of j such that 
p I (a + jqi.qi) and the corresponding value of M{q2^a-\- jq\) in the right hand side of (fTTI) 
would be undefined. The sum at the far right takes into account these primes. 

The validity of (fTOl) was checked immediately at the end of the computation of the constants 
M{q, a), for a fixed q and for every I <a<q with {q, a) = 1 by the same program that computed 
them. These results were collected in a file and a different program checked that (fTT)) holds 
within a very small error by building every possible relation of that kind for every qjE {3, . . . , 
100} and | <?2 with 1 <qi < qz- As in [7J, the total number of identities checked is 

100 100 

I I (p(^) = £(^-l-(p(^)) = 1907 

g=3 d\q q=3 
\<d<q 

but they are not independent on one another. We did not bother to eliminate redundancies since 
the total time requested for this part of the computation is absolutely negligible. Again as in 
[|71, the number of independent identities is 




= 1383, 



p<q 



where p denotes a prime in the sum on the left. Please remark that in [[Tl, page 323, we 
erroneously wrote that the previous sum is equal to 1408 which is in fact its value starting from 
n = \. 

Similar checks were done also for the 20 digits case. Working for every q < 300 we have 
12343 independent relations over a total number of 17453 ones. In this case, too, we obtained 
the desired precision (at least 20 decimal digits). 
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q 


a 




digits 


3 


1 


-0.3568904795094431291196495672231858954785 . . . 


104 


3 


2 


0.285054359023752579541743072498548421 1968 . . . 


104 


4 


1 


-0.2867420562261751986539451414394238573642 . . . 


104 


4 


3 


0.0482392690738179824093719800481197164157 . . . 


104 


5 


1 


-0.2088344499872589831393679436740355309848 . . . 


104 


5 


2 


0.3960964763519752181620428282992694487673 . . . 


104 


5 


3 


0.1386504040417767598465036287614177642579... 


104 


5 


4 


-0.26441521755885021 1 1 137516747779558229888 . . . 


104 


9 


1 


-0.1623582321428699054929449337179721787641 . . . 


104 


9 


2 


0.4073663127461732280783211701614365152217... 


104 


9 


4 


-0.1293374149143960665485300101 130823600639 . . . 


104 


9 


5 


0.0358267016686470538569841 87357 1 83 1423790 . . . 


104 


9 


7 


-0.0651948324521771570781746233921313566504... 


104 


9 


8 


-0.1581386553910677023935622850200712364040. . . 


104 


15 


1 


-0.1506479789635675321223227319951881293922. . . 


104 


15 


2 


0.396770207983 1 6025 19989220502990044 1 20830 . . . 


104 


15 


4 


-0.1298987796705018718274645424839556504587 . . . 


104 


15 


7 


-0.0006737316311850338368792219997349633156... 


104 


15 


8 


-0.1190129400473678821538466338276084167633 . . . 


104 


15 


11 


-0.058 1 8647 1 02369 1 45 1 1 70452 1 1 67884740 1 5926 . . . 


104 


15 


13 


-0.0756699892441886913329830707443071523119... 


104 


15 


14 


-0.1345164378883483392862871322940001725301... 


104 


21 


1 


-0.1084483613299595805404935908928381422038 . . . 


104 


21 


2 


0.4250487959922326653260015663325925353478 . . . 


104 


21 


4 


-0.1 122733018685413863141062477981428831803 . . . 


104 


21 


5 


0.1038169332452743207625287126777970125078 . . . 


104 


21 


8 


-0.0786267812146117135562454190553963721560... 


104 


21 


10 


-0.0827607926062097370238031413241567050640 . . . 


104 


21 


11 


-0.0174063652116240128169859448210067632915... 


104 


21 


13 


-0.0396915037660255136713627233956554116480... 


104 


21 


16 


-0.0990310666212991177281328116294234061275 . . . 


104 


21 


17 


-0.0495316505530113215988979410810678190789... 


104 


21 


19 


—0.0575425961745506509846081950401122043977 . . . 


104 


21 


20 


-0.0982465732345073585746579015543701721321 . . . 


104 



Table 1 . Some numerical results: the first column contains the modulus q, 
the second the residue class a, the third the computed value of M{q,a) and the 
fourth is the number of correct decimal digits we obtained. The table shows the 
values truncated to 40 decimal digits. 
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q 


a 




digits 


39 


1 


-0.0544150300747313383827161426970038945340 . . . 


104 


39 


2 


0.4598676271292146454635405190244096879502 . . . 


104 


39 


4 


-0.0459945989750685459192257387089186431466... 


104 


39 


5 


0.1419809783012313146832767050494466413940... 


104 


39 


7 


0.0795104580446772217182944213478570283508 . . . 


104 


39 


8 


-0.04827 1 127 16953638927878089938299013 1 8045 . . . 


104 


39 


10 


-0.0625407087664201913542100549743954739736 . . . 


104 


39 


11 


0.0351003379054991567928449805604752789995 . . . 


104 


39 


14 


-0.0462707514061093124433385415480332066056 . . . 


104 


39 


16 


-0.0671621258405927639177835818713861803563 . . . 


104 


39 


17 


-0.0045382608754434839604448757618956856127 . . . 


104 


39 


19 


-0.0078618584432586956665912002217299778289 . . . 


104 


39 


20 


-0.0459231434298270830511730272690723063199... 


104 


39 


22 


-0.0510428073342697409080043003007927978280 . . . 


104 


39 


23 


-0.01246316905345063503013701 13904825945643 . . . 


104 


39 


25 


-0.0581243810207459383640563019555017509681 . . . 


104 


39 


28 


-0.0517462190695606135794914158573956691966 . . . 


104 


39 


29 


-0.0270586216004553688028443336640559574303 . . . 


104 


39 


31 


-0.0298803394337857678868833935975175701503 . . . 


104 


39 


32 


-0.0473004404088889108385362179544297078717 . . . 


104 


39 


34 


-0.0463 158869206440559788 1 1065 1402015355032 . . . 


104 


39 


35 


-0.0548045111222561868022129723349733959006. . . 


104 


39 


37 


—0.0382400585981 196219570938701692763534202 . . . 


104 


39 


38 


-0.0652645592462251671904512528298502010370 . . . 


104 



Table 2. Some numerical results: the first column contains the modulus q, 
the second the residue class a, the third the computed value of M{q,a) and the 
fourth is the number of correct decimal digits we obtained. The table shows the 
values truncated to 40 decimal digits. 
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q 


a 


M{q,a) 


digits 


84 


1 


-0 


0734639142617973328342764883795225181917 . . . 


104 


84 


5 





1483235495915302335618054737355898008922 . . . 


104 


84 


11 





0290724024926249302145251081589204300848 . . . 


104 


84 


13 





0224290728747548696540000597026220130239 . . . 


104 


84 


17 





000300681 1962294747858026144190460423904 . . . 


104 


84 


19 


-0 


0057630788020940875837442254844273481232 . . . 


104 


84 


23 


-0 


0138497965201530520881 105236922199356889 . . . 


104 


84 


25 


-0 


06070167564290216083967764770082145 14160 . . . 


104 


84 


29 


-0 


0222266883543388763218294385273009788572 . . . 


104 


84 


31 


-0 


0325213918938471846809848382364695447795 . . . 


104 


84 


37 


-0 


0442703289531425263413021626262327744388 . . . 


104 


84 


41 


-0 


0475336593998265389573863164643172518681 . . . 


104 


84 


43 


-0 


034984447068 162247706217 1025 133 156240120 . . . 


104 


84 


47 


-0 


0445066163462559127992767610577927883843 . . . 


104 


84 


53 


-0 


04647876770424894303 1 5 1 1052979927 1 933763 . . . 


104 


84 


55 


-0 


0621205766407803833253627830982774246719 . . . 


104 


84 


59 


-0 


0498323317492407963847005555001138614693 . . . 


104 


84 


61 


-0 


0517795173724565634008639695556848562744 . . . 


104 


84 


65 


-0 


061 1014074876142825858879099751875289632 . . . 


104 


84 


67 


-0 


0515716262256392254744286000973214317642 . . . 


104 


84 


71 


-0 


0564000928602728372344159805280953932988 . . . 


104 


84 


73 


-0 


0502394007123625523428183030876871602845 . . . 


104 


84 


79 


-0 


0547607376681565913868306490031906316886 . . . 


104 


84 


83 


-0 


0507129138346808196172715850900529202640 . . . 


104 



Table 3. Some numerical results: the first column contains the modulus q, 
the second the residue class a, the third the computed value of M{q, a) and the 
fourth is the number of correct decimal digits we obtained. The table shows the 
values truncated to 40 decimal digits. 
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q 


a 




digits 


3 


1 


-0 


0179374320543395898017537423354360793084 . . . 


103 


3 


2 


-0 


2256492452247194384046517270072546894435 . . . 


103 


4 


1 


-0 


0303152628374217668471785632748622368557... 


103 


4 


3 


-0 


0922560086565230005866745667406677670593 . . . 


103 


5 


1 


-0 


0056989812258217866230186764730771910864 . . . 


103 


5 


2 


-0 


2072541594806942995597739906452261268831 . . . 


103 


5 


3 


-0 


0770818781394248684981698458665942749449 . . . 


103 


5 


4 


— 


002539881 893739366403827648178974475701 5 


103 


9 


1 


-0 


0020696391618847864572238027206860724807 . . . 


103 


9 


2 


-0 


1986304651091420386235919033135853530442 . . . 


103 


9 


4 


-0 


0039355992675986157162225504954464900241 . . . 


103 


9 


5 


-0 


0247398868156813518681399775437689383104... 


102 


9 


7 


-0 


01 19321936248561876283073891 193035168035 . . . 


102 


9 


8 


— 


00977889399989604791 991 98461 499003980889 

xjxj^z^ 1 ooy u ^jyyoy\j\j'~T i y L^y i yo'~T\jL'-Tyy\j\juyo\JOoy • . . 


103 


15 


1 


-0 


0007572379320997903470262134621931512735 . . . 


103 


15 


2 


-0 


1953264208891238586234409171756188321671 . . . 


103 


15 


4 


-0 


0016365552581033232050519597636557571823 . . . 


103 


15 


7 


-0 


01 19277385915704409363330734696072947160 . . . 


102 


15 


8 


-0 


0013342030920277845401475680955985955697 . . . 


103 


15 


11 


-0 


0049417432937219962759924630108840398129 . . . 


102 


15 


13 


-0 


0036159002725660353133424956399798761365 . . . 


102 


15 


14 


— 




103 


21 


1 


-0 


0003412956292374148069148220346920460252 . . . 


103 


21 


2 


-0 


19423449473346996888940039745601 12287700 . . . 


103 


21 


4 


-0 


0002098816767539160024207752141222498501 . . . 


103 


21 


5 


-0 


023509322852284191 1201270554154748282363 . . . 


102 


21 


8 


-0 


0007915893971472685099562881470728940934. . . 


102 


21 


10 


-0 


0006922252022738492137317334872558725556 . . . 


102 


21 


11 


-0 


0046535617744727410631965497667779194062 . . . 


102 


21 


13 


-0 


003248 13 1763552583 193 1756409952984059420 . . . 


102 


21 


16 


-0 


0004945238101126788996808520807420758727 . . . 


102 


21 


17 


-0 


0020363536140986162458739759280171063102 . . . 


102 


21 


19 


-0 


0016578370022937005358116763179908292912. . . 


102 


21 


20 


-0 


0004239228532466525760974602939007126272 . . . 


102 



Table 4. Some numerical results: the first column contains the modulus q, 
the second the residue class a, the third the computed value of B{q,a) and the 
fourth is the number of correct decimal digits we obtained. The table shows the 
values truncated to 40 decimal digits. 



10 



q 


a 




digits 


39 


1 


-0.0001 121391210993880688819721271925997627 . . . 


102 


39 


2 


-0. 193476965597537 1 9934908 1 37692 1 06 1 9094240 . . . 


102 


39 


4 


-0.0002995815105830464353463369791216808136 . . . 


102 


39 


5 


-0.0232346770448918237834268807192861283189 . . . 


102 


39 


7 


-0.01 13299371 144385520676547127480047428604 . . . 


102 


39 


8 


-0.0002454341727558259212509705166966869744. . . 


102 


39 


10 


-0.0000459138224155805565513538531922246471 . . . 


102 


39 


11 


-0.0044929480570028872659134897718135000937 . . . 


102 


39 


14 


-0.0002193611514478558639221844642199225052. . . 


102 


39 


16 


-0.0000223235704255911196145171142624220894. . . 


102 


39 


17 


-0.0018338873065642236614942324023402025533 . . . 


102 


39 


19 


-0.0015031092616855056695375513171955283284... 


102 


39 


20 


-0.0001858990560853595070479345313080934723 . . . 


102 


39 


22 


-0.0001705919794542183214877183738499262132 . . . 


102 


39 


23 


-0.00 10537 1 054582 1 54353243 8097 1 90 1 947683990 .. . 


102 


39 


25 


-0.0000730986515003499343341780315856293861 . . . 


102 


39 


28 


-0.0001342213885419137800239278870137043470 . . . 


102 


39 


29 


-0.0006675737489094563083519829970440655189 . . . 


102 


39 


31 


-0.00058526763 15 15 17 17436339529 174197209276 . . . 


102 


39 


32 


-0.0001433829254050606314390524017643892945 . . . 


102 


39 


34 


-0.0001404190990760282183009755731128149793 . . . 


102 


39 


35 


-0.000070729147470948299198995013088791 1099 . . . 


102 


39 


37 


—0.0004011981531447411395316606101222440231 . . . 


102 


39 


38 


-0.0000246764708272542810865300784362317789 . . . 


102 



Table 5. Some numerical results: the first column contains the modulus q, 
the second the residue class a, the third the computed value of B{q, a) and the 
fourth is the number of correct decimal digits we obtained. The table shows the 
values truncated to 40 decimal digits. 
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q 


a 


B{q,a) 


digits 


84 


1 


-0.0000119163858637686167954725330316682793 . . . 


102 


84 


5 


-0.0232403602184713008048627754438543014690... 


102 


84 


11 


-0.0044365093956013183002165422530420512818 . . . 


102 


84 


13 


-0.0032012002462998617358456975415447057261 . . . 


102 


84 


17 


-0.0018681454567877949532487758242611745268 . . . 


102 


84 


19 


-0.0014999868517105941255280000567662988513 . . . 


102 


84 


23 


-0.0010419325887589517242899271 1 1 1822879451 . . . 


102 


84 


25 


-0.0000667248398877511948366874331992186367 . . . 


102 


84 


29 


-0.0006749316281200272421803236576904643492 . . . 


102 


84 


31 


-0.0005598064648252633420494915697566166895 . . . 


102 


84 


37 


-0.0003842227265021586749603313599748870266. . . 


102 


84 


41 


-0.0003 1 505848382028083 1 624224902 1 937896620 . . . 


102 


84 


43 


-0.0003293792433736461901 193495016603777458 . . . 


102 


84 


47 


-0.0002689626338128903152642799716205267673 . . . 


102 


84 


53 


-0.0002170523788714227629800075137358681243 . . . 


102 


84 


55 


-0.0000469315172527214573299434537537002158 . . . 


102 


84 


59 


-0.0001682081573108212926252001037559317834... 


102 


84 


61 


-0.0001578501505831064102836762612245304399 . . . 


102 


84 


65 


-0.0000453815847657077478783488866523727492 . . . 


102 


84 


67 


-0.0001431568368661648075840877809230312134... 


102 


84 


71 


-0.0001 166577690272412677759644893824297442 . . . 


102 


84 


73 


-0.0001324187374485858716822419174992558660. . . 


102 


84 


79 


—0.0001103010836105202247205207207671888461 . . . 


102 


84 


83 


-0.0001088643694263717444732353917069229651 . . . 


102 



Table 6. Some numerical results: the first column contains the modulus q, 
the second the residue class a, the third the computed value of B{q,a) and the 
fourth is the number of correct decimal digits we obtained. The table shows the 
values truncated to 40 decimal digits. 
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